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Our member from the first days, Herbert Appel asks for publishing 
the following announcement. I hope I can help you and much luck, 
Herbert . 


NEW - MATH - VERLAG 

Der Verlag fur den innovativen Mathematikunterricht 


Wir bieten: 

Nach bewahrten didaktischen und methodischen Gesichtspunkten ausgearbeitete Komplettlosungen zu 
Themenbereichen der Sekundarstufe I und II fur den Einsatz der Programme CABRI GEOMETRE 
(Zug-Geometrie) und DERIVE in Ihrem Mathematikunterricht. 

Wir suchen: 

Autoren, die iiber Unterrichtsmaterialien zum Einsatz o.g. Softwareprodukte verfiigen. Geeignete Lo- 
sungen werden wir nach eingehender Priifung in unser Angebot aufnehmen. 

Informieren Sie sich! 

New-Math- Verlag, Am Weiher 5, D-8643 Klips, Tel: 09264-8853, Fax: 09264-8799 


Some food for your memory: DERIVE XM 

DERIVE XM operates identically to the regular DERIVE but can take 
advantage of up to 4 gigabytes (4 billion bytes) of extended memory. 
The size of problems that DERIVE XM can handle is now limited only 
by your computers ' s memory and speed. 

Operations such as expanding polynomials, inverting matrices, and 
solving sysems of linear equations can require a great deal of mem- 
ory to store the intermediate and/or final results. For example, the 
maximum size random number matrix that regular DERIVE can invert is 
about 40 by 40. However, DERIVE XM was able to invert a 100 by 100 
matrix an a computer with 4 MB of total memory (i.e. of extended me- 
mory) . 

However, the system requirements for DERIVE XM are still rela- 
tively modest: a 386 or 486 based PC compatible computer with a mi- 
nimum of 2 MB of memory (i.e. 1 MB of conventional and 1 MB of ex- 
tended memory) . 


Use MathPORT to Drive DERIVE to the Publisher: 

With DERIVE and MathPORT you not only get the 
right answer, but you can share your professional- 
quality documents with others! MathPORT is the 
leader in converting math equations for word proc- 
essed documents. 

The equations beneath were created and ublished as 
TIFF graphics files by MathPORT. They were auto- 
matically converted from DERIVE files without 
editing. 

(K-Talk Communications 
30 West First Ave . Suite 100 
Columbus, OH 43201) 
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Liebe DERIVE Anwender! 

Ich mochte Sie alle zu Beginn des 3. DNL-Jahres 
recht herzlich begrul.Sen und Ihnen auf diesem Weg 
fiir die vielen lieben Grii lie anlasslich des Jahres- 
wechsels danken. Die breite Zustimmung fiir die 
Ge-staltung des DNL hat mich sehr gefreut. Viele 
von Ihnen haben den kleinen Fragebogen aus 
DNL#8 ausgefiillt und so ist ein ganzer Berg von 
Anregungen zusammengekommen. Da noch immer 
Rtickmeldungen einlangen, mochte ich erst im 
nachsten DNL eine Zusammenstellung Ihrer Ideen 
und Anregungen veroffentlichen. 

Natiirlich mochten einige von Ihnen mehr schulbe- 
zogene Themen, andere befiirchten eine allzu grol.Sc 
Schullastigkeit. Einige sind Grafikliebhaber - so wie 
ich, ich kann das nicht verleugnen - andere empfin- 
den manche Dinge als Spielerei. So werde ich nun 
versuchen den Wiinschen halbwegs gerecht zu wer- 
den. 

Eine besonders nette Geschichte hat der Artikel 
von Marko Horbatsch: Ein Grazer DUG Mitglied 
schickte mir eine Kopie, darauf hin nahm ich in 
meinem „besten Englisch“ Kontakt mit M. Horbatsch 

auf und bekam sehr rasch einen netten Brief in 

bestem Deutsch zuriick. Herzlichen Dank Herr Hor- 
batsch fiir Ihre Erlaubnis zur Veroffentlichung. 

Fiir den nachsten DNL werde ich einem haufig 
genannten Wunsch nachkommen und den Schwer- 
punkt Statistik und Wahrscheinlichkeitsrechnung 
setzen. Vielleicht hat jemand von Ihnen einen Bei- 
trag zu diesem Thema auf Lager? 

Vielfach wurde der Wunsch nach einer Informati- 
on iiber die Moglichkeiten der neuesten DERIVE 
Versionen geauGert. Sie konnen auf Seite 41 eine 
Zusammenstellung der wichtigsten neuen Features 
finden. (2.50 und hoher) Beachten Sie auch den 
kurzen Bericht in den benachbarten Informationen 
liber DERIVE XM. 

Aufgrund der - an sich erfreulichen - Materialfiille 
muss ich Sie bitten, fiir dieses Mai auf den BBS zu 
verzichten. Auch das User Forum wurde sehr kurz 
gehalten. Einige umfangreiche Briefe warten auf 
ihre Veroffentlichung. 

Absch I iel.Scnd bitte ich alle jene, die ihren Mit- 
gliedsbeitrag fiir 1993 noch nicht entrichtet haben, 
dies bei nachster Gelegenheit nachzuholen. Herzli- 
chen Dank. 

Bis zum nachsten Mai mit den besten Grii l.icn Ihr 


Dear DERIVE Users, 

First I want to welcome you to the start of DUG's 
third year of existence. I also want to thank ever- 
ybody for the nice greetings and for the many good 
wishes I and my wife received for the New Year. 
Many of you have sent back the little form from 
DNL#8 and so a lot of ideas have assembled in my 
DNL-folder. So we are receiving answers every day 
I want to give a list of the ideas, suggestions and cri- 
ticisms in the next DNL. 

Of course some of you want to have more con- 
tributions connected with didactics and teaching 
mathematics while others may be afraid of too much 
schoolmathematics, they would prefer more scienti- 
fic applications. Some like graphics - as I do, I can't 
deny it - others dislike them as a triviality. I'm trying 
to find a way to satisfy most of your interests 

Marko Horbatsch's "Physics in the Classroom" has 
an especially nice story. A former DUG member 
from Graz, A, sent me a copy and I contacted 
M. Horbatsch in Canada in my 'best English' and 

very soon after I received a friendly letter 

but in "best German" ! Many thanks Marko for gi- 
ving your permission to reprint your interesting 
work. 

In the next DNL I will fulfill a wish mentioned 
from many members and concentrate on statistics 
and probability theory. Maybe there are some read- 
ers who have a DERIVE-application on these sub- 
jects? 

Widely I've noted the wish for information about 
the facilities of DERIVE’s latest versions. On page 
41 you will find a list of the most important new 
features (2.50 and higher). Please note the comment 
on DERIVEyXM on the ’Information page’. 

Because of the - really delightful - quantity of 
material received I must ask you to do without the 
BBS this issue. The User Forum had to be short- 
ened, too. Some extensive letters are waiting for 
publication. 

Finally I’d like to ask the DUG-Members who 
haven't paid their membership fees for 1993 to do 
this at their earliest opportunity. Thank you. 

Until next time 
with my best regards 





Contributions : 


The Derive-News-Letter is the 
Bulletin of the Derive-User 
Group. It is published at least 
three times a year with a con- 
tents of 30 pages minimum. The 
goals of the D-N-L are to en- 
able the exchange of experi- 
ences made with Derive as well 
as to create a group to discuss 
the possibilities of new me- 
thodical and didactic manners 
in teaching Mathematics. 

Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/ FAX: 43- (0) 2275/8207 
e-mail: no j o . boehm0pgv . at 


Please send all contributions 
to the above address . Non- 
English speakers are encouraged 
to write their contributions in 
English to reinforce the inter- 
national touch of the D-N-L. It 
must be said, though, that non- 
English articles will be warmly 
welcomed nonethe-less . Your 
contributions will be edited 
but not assessed. By submitting 
articles the author gives his 
consent for reprinting it in 
the D-N-L. The more contribu- 
tions you will send to the Edi- 
tor, the more lively and richer 
in contents the Derive-News- 
Letter will be. 


Preview: (Contributions for the next issues) 

Census Assigment, Some Worksheets, K. Eames, UK 

Bisection with DERIVE , D.M. Dyer, USA 

Fluid Flow in DERIVE , Reuther a.o., BRA 

Newton's Chaos, Graphic Integration, J. Bohm, AUT 

Computer Aided Mathematics in School, K.H. Keunecke, GER 

Continued Fractions, R. Setif, FRA 

DERIVE in Hawaiian Classrooms, Sawada, USA 

Stability of ODEs, Kozubik, SFK 

Bivariate Normal Distribution, Marinell, AUT 


DNL#10 will be published June 93 


Download all DNL-Derive- and Tl-files from 

http : / / www . austromath . ac . at/ dug/ 
http : // www . bk- teachware . com/ main . asp?session=37 505 9 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D’Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag Josef Bohm 
Herstellung: Selbstverlag 
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Reuther & Domingos, Petropolis, Brazil 

We would like to submit our preprint "Extended Fourier Series with DERIVE" to the DNL. We hope 
it will fulfill all the specifications. 

We would also like to remember that this is our second participation, because a member of our group, 
Patricia R. Medici, has been in the congress happened in Krems in April. Our group is a small one, 
and we're trying to widespread the DERIVE "culture" in Brazil. 

K.Eames, Kingsway College, UK 

On receiving the recent copy of DNL I thought you might like to know how DERIVE is being used at 
Kingsway College, which is a Further Education College. The reasons why DERIVE is being looked 
at are as follows: 

1 . Kingsway is interested in the impact of such packages on the Mathematics curriculum, I felt 
DERIVE is better for the 16-19 age group than REDUCE or MATHEMATIC A, the former seemed to 
be too complicated and the latter too memory hungry. 

2. Kingsway is already producing manuals and assignments for A-level and University of North Lon- 
don students to try out on their respective courses. Student manuals on DERIVE's Basic commands 
and Use of DERIVE in Algebra complete with simple test exercises for students to complete, in order 
for them to see if they have mastered the relevant commands. There are also manuals on DERIVE's 
graphical capabilities along with worksheets and on DERIVE's use in Calculus. 

3. Kingsway is currently looking at the possible use of DERIVE as a means of introducing certain to- 
pics, Matrices for example. 

plus I have given a number of INSET sessions on the Use of DERIVE for the London Branch of the 
Mathematical Association and for Suffolk College of Higher Education. A copy of the INSET manual 
is enclosed. Many of the sheets arc being updated at this present moment in time. 

My own views on the use of DERIVE arc as follows: 

1 . In supporting investigative work, thereby enabling more time to be spent on formulating problems 
and predecting possible results 

2. To introduce topics 

3. To enhance students opportunities in their course and to come into contact with it before going onto 
HE 

4. To check results to problems solved without the use of DERIVE. 

K.de Ridder, Zwolle, NL 

Wij zijn met een experiment bezig om computeralgebra te introduceren oop de H.T.S. Zwolle, studie- 
richting elektrotechniek. Dat experiment begint februari 1993. Het gebruik van CA als wiskundig ge- 
redschap vor het oplossen van technische problemen staat darbij voorop.(the use of CA as a mathema- 
tical tool to solve technical problems is the most important aim). Als U geinteresseerd bent willen we 
U graag op de hoogte houden van onze vorderingen en resultaten. 


DNL: I think we would be interested in your experiments and we hope to hear from you when 
you have got any results. 

There are some other extensive letters - from MryNeurath, Mr Propper, Mr Schmidt, Mr Wil- 
liams, Mr Horbatsch and others - which I want to publish in the next issue. Many thanks for 
your interest and cooperation. 
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Solving Physics Problems in the Classroom with DERIVE 2.0 

Marko Horbatsch, York, Ontario, Canada 


Abstract: We present an overview of problems that can be solved now in DERIVE 2.0 due 
to the incorporation of routines that handle iteration as well as other programming elements 
such as logical branching. The applications include recursive calculation of orthogonal po- 
lynomials, solution of iterative maps, numerical solution of ordinary and simple partial dif- 
ferential equations. 

DERIVE ® A Mathematical Assistant [1] has proven to be a useful educational as well as research 
tool [2-4]. With the new version 2.0 some of the major constructive criticisms raised in our previous 
review have been taken into account. The basic philosophy of the program reamins the same, namely 
to achieve compact code that can be loaded on almost any MS-DOS based machine irrespective of the 
CPU. This is achieved by defining intrinsically to the program only a minimal set of functions and 
relegating all special functions to external files. The improvements to the package are in two areas. On 
the program side one has added to the symbolic environment two important aspects of a programming 
language, namely conditional branching as well as iteration. On the support files new files with func- 
tion definitions were added and a description of these is provided in a much enhanced manual [1], 
Also the intrinsic help utility covers now the utility files in a brief fashion so that it is possible to use 
them without the manual. Extensions to the already ample graphing capabilities of the program include 
the possibility to draw user calculated isometric projections of three-dimensional curves, surfaces and 
vectors in two-dimensional plot windows as well as graphing of complex functions. 

The putpose of this review is to demonstrate how Derive 's new features can be combined with the 
graphical capabilities to extend significantly the program's usefulness for educational and research 
putposes. It will be shown using such examples as dicrete maps, solution of ordinary and simple par- 
tial differential equations that it is possible to set up a theoretical physics laboratory with this envi- 
ronment. 

A major problem for the instructor of a course in modern physics in an undergraduate curriculum 
is to familiarize the student with the idea of a differential equation - eigenvalue problem. Instead of 
concentrating on the physics behind Schrodinger ' s equation the student is preoccupied with learning 
the mathematics. The symbolic -graphic environment that Derive provides is ideally suited to present 
results and allow self-study of differential equations, orthogonal polynomials, their recursion relations, 
etc. 

Special software packages are in existence for putpose of plotting potentials and the eigenvalues 
and eigenfunctions of simple Hamiltonians. The following example will show ow the same can be 
achieved with a few programming steps in Derive. An advantage over dedicated packages is that the 
student has to understand how Derive managed to obtain these results. 

In detail: Line #1 in fig. la provides a definition for the Hermite polynomials, which is included 
in one of the support files that come with the software. Derive provides access to many special func- 
tions in this way, i.e. one can use after importing such a line of code into memory. While this may 
seem awkward compared to programming languages as MAPLE, where such functions can be called 
as intrinsic functions, for educational putposes, one has the advantage of seeing a definition of the 
function displayed. In the new version of Derive utility files can also be loaded into memory without 
displaying them as a part of the user's program. 
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# 1 : 


n x 


HERMITE_H(n, x) (-1) -e \dxj 


n 2 

- x 

e 


PSICx, n) :: EXP(- 0.5-x ) • HERMITE_H(n, x) ■ 


*2: 


Jtt- 2 • n ! 


#3: 

HAMILTON(u) := 0.5- 

r 


2 2 'l 

U + X -u 




^dxj 

J 


#4: 


ev 


r 

VECTOR ■ 


HAMILTON(PSICx, 
PSICx, i) 



i, 0, 3 

J 



#6: VECTORCO. 5 • PSI(x, k - 1) + ELEMENT(ev, k), k, 1, 4) 


Figure la: A calculation of the lowest four eigenvalues and eigenfunctions of the 
harmonic oscillator (in units h-m-w- 1). In line #4 DERIVE calculates the ei- 
genvalues with results shown in line #5. Line 6 provides a vector in which the low- 
est four eigenfunctions shifted by the corresponding eigenvalues are generated . 

Line #2 defines the n h harmonic oscillator eigenfunction with appropriate normalization factor. 
As Derive has no intrinsic knowledge of the Hermite polynomials and their properties, it is not able to 
prove in general that the (real) functions defined in this line satisfy the orthonormality condition 

I V' n (x)yr m {x)dx = S nm . (1) 

J —00 

It is, however, possible to verify this property for individual cases of m and n and this is calculated 
exactly in Derive. 

Line #3 provides the Hamiltonian as it acts on some arbitrary function u, which for our choice of 
potential and units reads: 


H 


Id I 2 1 2 

- + -x- 

2 dx 2 2 


( 2 ) 


Note that the operation with the Hamiltonian operator on u is defined as a function with u as an argu- 
ment., but that this function includes the action of a differential operator on u. Line #4 calculates a 
vector that contains the lowest four eigenenergies form the equation Hy/= Ey/, with the result dis- 
played in line #5. Lor plotting purposes a vector is generated whose elements contain the lowest four 
eigenfunctions shifted from each other by the eigenvalues. This is done in order to show in fig. lb the 
spatial extent in which a classical one-dimensional harmonic oscillator would move, if it was given the 
corresponding amount of energy. 




Figure lb: The eigenfunctions displayed in a Derive 2D-Plot Window 


By taking the square of the eigenfunctions one could show on a similar graph the probability of find- 
ing a quantum mechanical particle at a given location x (assuming that it is in a given eigenstate n) and 
draw attention to the oscillatory structure of this quantity as well as the classically forbidden region. 

One could argue that such a figure is included in many textbooks on modern physics or quantum 
mechanics and ask whether a student will learn much from seeing the same on a computer screen. 
Furthermore, the essentially effortless calculation summarized by four lines of Derive code could dis- 
courage students from performing these calculations by hand. An appropriate way to make use of De- 
rive for students who are supposed to master the mathematics involved, would be in such a way that 
the program would be used for checking and gaining confidence in a calculation by hand. Solutions to 
radial Sclirodinger equations for several commonly used potentials involving other orthogonal poly- 
nomials can be obtained in the same way. 

The calculations shown so far were feasible in previous versions of Derive. A novel feature is the 
possibility to study recursion relations for symbolic expressions, one of the strong points of symbolic 
languages. The line of Derive code given below provides an alternative calculation of the Hermite 
polynomials based upon the recursion relation 

H m+1 (x) = 2xH m (x)-2MH m _ 1 (x), (3) 

which is started by H 0 (x) = 1 and H\(x) = 2x. 

The new function in Derive 2.0, which generateds recursions is ITERATE. Its syntax is given by 
ITERATE (/(w), w, w 0 , n) and it returns the n lh iteration of the assignment w —> ft w) starting with w 0 - A 
similarly working function ITERATES returns all elements of the iteration. The first three arguments 
are allowed to be vectors, which in our case enables the two-step recursion to be calculated. In our 
example the symbolic placeholder w is a three-vector as we have to update the iteration index in addi- 
tion to storing H M (x) and H M - i(x). Reference to the individual elements in the recursion is made by 
Derive 's routine to extract elements from vectors (and matrices), ELEMENT (v,j) returns the f h ele- 
ment of vector v. Thus the recursion in (1) is programmed in Derive by the following line to generate 
H n (x). 
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ELEMENT (ITERATE ( [ELEMENT (w, 2), 2 *x*ELEMENT (w, 2) - 

2* (ELEMENT (w, 3) + 1 ) *ELEMENT (w, 1), ELEMENT (w, 3) +1], w, 
[1, 2*x, 0], n) , 1) 


Here the first element of w returns H n (x), the second stores //„_ i(jc) and the third element of w contains 
the index n. The square brackets in Derive code refer to a vector given in terms of its elements. The 
fact that one has to update the iteration index and refer to it by the ELEMENT -function is clumsy and 
this should be simplified in some future release. Once the student has learned how to set up recursions 
in Derive, he can save a lot of time by avoiding tedious calculations, compare results with the defini- 
tion given in line #1 of fig. la and concentrate on the content rather than the technical aspects of the 
calculation. 


(ITERATE (Tw , 2-x-w - 2-(w + l)-w , w + ll, w, [1, 2-x, 0], rQ) 

L 2 2 3 1 3 J 

1 

Recent Derive versions don't need the 'clumsy' ELEMENT-function 

This new capability of handling iteration opens up possibilities held traditionally by numeri- 
cal/graphical environments, such as Quickbasic or Turbo-Pascal. Derive is easier to program than 
these environments and has an easier access to graphing in multiple windows. We have chosen an 
example from classical mechanics that cannot be solved in closed form. Its intent is to illustrate possi- 
ble trajectories for a classical charged particle with small mass that impinges on a two dimensional 
fixed array of oppositely charged ions. 

The ITERATES function is used to integrate Newton ' s equation for a particle moving through an 
array of point-charges. The dimensionless equation (m = e = 1) is solved approximately by the so- 
called leapfrog algorithm coded in line #11 below: 

y (u) = v (U 2 ) + ^F(r(f„)), r{t n+l ) = r{t n ) + dt\{t n+U2 ). (4) 


Here r, v and F are the position, velocity and force vectors respectively and r(t 0 ) and \(t \a) the initial 
conditions. This fixed time-step algorithm is not very accurate but economical for many applications. 
In order to avoid accuracy problems, which would force one to use a very small time step dt, the po- 
tential of the static sources has been regulated: while the motion of the test charge is confined to the x- 
y plane, the external sources can be thought of as resinding in a plane located at a distance x = yfa 
above (or below). The solution generated by LEAP (n) is a vector variable named vek, which con- 
tains the *- and y-coordinatcs of r(f 0 ) and One can see that the coding is somewhat clumsy due 

to the fact that one has to iterate on a single object and refer to its elements. The code shown below 
would appear less intimidating, if we had changed the ELEMENT function to one with a shorter name. 

The routine EXTRACT_2_COLUMNS is supplied in one of the support files and serves to ex- 
tract two colurns out of a given matrix, which is achieved by transposing the matrix, extracting two 
rows and transposing again (lines #1 and #2). xvek and yvek are vectors denoting the x- and 
y-coordinatcs of the points at which the fixed charges are located ncent calculates the number of these 
charges, areg is used to regulate the interaction and lines #7 and #8 provide the x- and y-components 
of the force acting on a test-charge located at (x_, y_) . incond specifies the initial condition 
r 0 = (2,0.4) and Vo = (-0.5,0) of the test-charge in the plane. As the ITERATES function returns a 
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matrix containing positions and velocities on the time mesh one has to the columns containing the 
(x,y)- values in order to graph the trajectory. 


#1: 

EXTRACT_2_ELEMENTSCv, p, 

q) := [ELEMENTCv, p) , ELEMENTS, q)] 

#2: 

EXTRACT_2_C0LUMNSCa, cl, 

c2) := EXTRACT_2_ELEMENTSCa' , cl, c2) ' 

#3: 

xvek := [-1, 0, 1, -1, 0, 

1, -1, 0, 1] 

#4: 

yvek := [0, 0, 0, 1, 1, 1, 

-1, -1, -1] 

#5: 

ncent := DIMENSIONCxvek) 


#6: 

areg := 0.5 



ncent x_ - ELEMENT(xvek , i) 

FX(x_, y_) := I 

#7: i=l 2 23 

7((x_ - ELEMENT(xvek, i)) + (y_ - ELEMENTCyvek, i)) + areg) 

ncent y_ - ELEMENTCyvek, i) 

FY(x_, y_) := I 

#8: i=l 2 23 

V((x_ - ELEMENT(xvek, i)) + (y_ - ELEMENTCyvek, i)) + areg) 

#9: incond := [2, -0.5, 0.4, 0] 

#10: dt := 0.01 

#11: LEAPCn) := ITERATESC1 ELEMENT Cvek , 1) + dt- ELEMENTCvek, 2), ELEMENTCvek, 2) - 

dt-FXCELEMENTCvek, 1), ELEMENTCvek, 3)), ELEMENTCvek, 3) + dt- ELEMENTCvek, 4), 
ELEMENTCvek, 4) - dt- FYCELEMENTCvek, 1), ELEMENTCvek, 3))], vek, incond, n) 

#12: EXTRACT_2_COLUMNSCLEAPC1200) , 1, 3) 



Figure 2: Numerical solution of Newton's, equation for a test charge impinging on 
an array of fixed charges located in a plane with the leapfrog algorithm. 

(calculated with DERIVE XM within 1347 sec; DERIVE 2.54 wouldn't do it in one run, I 
had to split. DERIVE 6.10 needs 7 sec in 2006. Editor) 


See the leapfrog algorithm - performed with DERIVE 6 - in the following. The output looks the same! 
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#1: pts 


-1 0 

0 0 

1 0 

-1 1 

0 1 

1 1 

-1 -1 

0 -1 


1 -1 


#2: [areg := 0.5, incond := [2, -0.5, 0.4, 0], dt := 0.01, vek =] 


*3: 


DIM(pts) 

FX_(x_, y_) := I 
irl 


x_ - pts 


i .1 


2 2 ^3 

(x_ - pts ) + (y_ - pts ) + areg 

i , 1 i ,2 > 


#4: 


#5: 


DIM(pts) 

FY_(x_, y_) := I 
i=l 


pts 


i ,2 


2 2 ^3 

(x_ - pts ) + (y_ - pts ) + areg 

i , 1 i ,2 > 


LEAP_(n) :r ITERATES ( 


[vek + 

L 1 


dt-vek , vek - dt-FX_(vek , vek ), vek + dt-vek , 


vek - dt-FY_(vek , vek )j, vek, incond, n) 


#6: 


(LEAP_(1200)) COL [1, 3] 


The results shown in fig. 2 are produced 
upon the simplification of line #12 and dis- 
played as a sequence of points. These points 
indicate the position of the moving charge at 
discrete time-steps. Subsequenent points are 
displayed in different colors chosen from a 
palette on the computer screen and could be 
represented in the printout by different grey 
shades. From the graph one obtains a feeling 
for the speed of the test charge in the various 
parts of the trajectory. The spacing of the 
points indicates where the probe charge accel- 
erated and decelerated, etc.. Further interesting 
quantities could be calculated and displayed as 
function of time. One can easily study the sen- 
sitivity of such a system on the initial condi- 
tions and watch how a minimal 


change in the initial position or velocity makes the 
test charge take a different turn at some point during 
its path. 
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(LEAP_(2400)) COL [1, 3] incond := [1.5, -0.3, 0.5, 0.2] 




A classical mechanics laboratory can be set up following this example, i.e., one can assemble a 
number of problems with different potentials of relevance to physics, chemistry, etc., for numerical 
study. An interesting comparison comes to mind with existing numerical software packages dedicated 
to such problems, such as Chaos Demonstrations (published by the American Physical Society), Chaos 
Dynamics Workbench (APS), Phaser (Springer) and others. On this front Derive cannot compete at all 
in terms of speed, as all computations are done without a coprocessor and usually in exact precision 
(rational) arithmetic. Another technical problem is the fact that all values of the iteration are stored and 
extracted for plotting only subsequently. If one is forced to a small time-step for the integration, this 
may lead to storage problems due to the DOS 640 kB memory barrier. However, from an educational 
point of view Derive is preferable, as, again, one is forced to understand what is going on. For the 
production of nice figures, where large amounts of data and plotting on the run are required, one could 
then turn to one of the programs mentioned above. 

At this point it should be clear that Derive provides an interesting tool for the numerical analyst 
(as long as one is content with basic problems, such as occur in the classroom). Without having to deal 
with round-off errors (if the problem is solvable in the rational numbers) one can study the effects of 
dicretization errors, and display results conveniently. However, it is practically impossible to write 
complicated programs, as Derive still cannot be called a programming language. (Things have 
changed, Ed.) While many elements of looping, conditionals, iteration arc now implemented, one has 
to live with the structure of single line calls. To code a complicated branching expression is possible, 
but such a single line can become quite long and very difficult to debug. (I repeat: Things have 
changed, Ed.) 

There exist examples for which the VECTOR function which covers one element of looping, will 
not simplify its main argument. I found it impossible, e.g., to generate a sequence whose elements 
were determined by the approximate solution of a nonlinear equation, even though individually De- 
rive's SOLVE function calculated each element very quickly. From the point of view of symbolic 
language syntax Derive 's relationship to a complete language (such as Maple, Reduce, etc.) can be 
compared to that of early Basic dialects to Fortran. On the other hand it is so much easier to use than 
any of the fullblown symbolic programming languages known to me that it should definitely be used 
as an introduction to symbolic computing. 

Marko Horbatsch, Dept of Physics & Applied and & Computational Mathematical Science 
Programme, York University, Toronto, Ontario, M3J 1P3, Canada 
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#18: VECTOR(herm(k) , k, 0, 10) - VECTOR(HERMITE_H(k) , k. 0, 10) 
#19: [0, 0, 0. 0, 0, 0. 0, 0, 0, 0, 0] 


e •henn(4)-hen«(6) dx = 0 


#21: - x 

e •herm(13)*herm(27) dx = 0 


#22: - x 

e .herni(4).herm(4) dx = 384 .,/n 


VECTOR J e •herm(k)-herm(k) dx, k, 0, 5 


#24: Ijn, 2 ‘Jtt, S-Jn, 48-^, 384 '-Jn, 3840-,/7t] 

k 

#25: VECT0R(2 - k ! - ^rr , k, 0, 5) = [Vtt, 2-Jn, S-Jn, 48-Vtt, 384. Jn, 3840 ■■Jn] 

It could be a nice challenge for the students to find out the general rule to generate sequence #24 


#26: - x 

e •lierm(15).henn(15) dx = 42849873690624000-> 


#27: 2 -15! •> = 42849873690624000 •> 

In addition to the recurrence equation given in Marko Horbatsch 's contribution the Hermite Polyno- 
mials follow another recurrence equation which is connected to a differential equation: 

H"(x) - 2 xH' n (x) + 2 nH(x) = 0 


ffd )2 d 

VECTOR — herm(k) - 2-x — herm(k) + 2*k»herm(k), k, 10, 15 
AdxJ dx 


[ 0 , 0 . 0 . 0 , 0 , 0 ] 


H' ll (x) = 2nH n _ l (x) 


VECTOR — herm(k), k, 10, 20 - VECT0R(2-k-herm(k - 1), k, 10, 20) 

Idx J 


[ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ] 


Hermite Polynomials are also defined using the weight function co(x) - e 2 . 
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Last July I received a letter from Finland 

"What do you like my files? The research in mathematics and phy- 
sics (mechanics) has been long time my work, although I am not a 
professional physicist. I am a writer (middle level studies in phys- 
ics, too), and I have new ideas about mechanics. Are the files suit- 
able for circulation? The mathematics in these files is very sim- 
ple . " 

I'm not a physicist, too, so I can't give a comment on Mr Ahonen's ideas. The time he has 
spent is as remarkable (a lot of hours for calculating and plotting) as his enthusiasm. When I 
wrote to Finland, that I would like to publish his contribution I received 5 versions of Ahonen's 
Mechanics, one time two letters a day. In December when I informed Mr Ahonen that his 
ideas would be published in March 93 he seemed to be a bit disappointed: "So much work 
and now so long to wait! This is for me the most important thing, and for the world, too". 

But finally, here is 

The Mechanics of Erkki Ahonen 

(Part 1) 

The Mechanics of Erkki Ahonen 

Main Features 

Erkki Ahonen Salo Finland 

[DERIVE V.1.22 used. V.2.51 does not work well] 


I have been working at the basic mechanics many years and have got 
interesting results. The fundamental result of my research is the 
equation of motion 



where a = d 2 x/dt 2 = du/dt = u du/dx = acceleration, 
a 0 = constant acceleration; 

u = dx/dt = velocity, u 0 = constant velocity; p,k,n e R 
d = real or imaginary number. 

x = coordinate of place, t = coordinate of time. 


The integration constants x 0 ,t 0 move the origin. 
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EXAMPLE: The equation of the parabola is 

a o( t - t o) = ± V 2 a o( x - x o) / 
when p,n = 0 in Eq. (1) : a = ao . 

Eq. (1) gives as solutions a large set of curves when the 
parameters are changed. 

DYNAMICS: The equation of force is 

f(t) = m 0 , ( 

v ' dt bL 


where m 0 = rest mass; bL=^l-|^— J ;; 
c = velocity of light. 


With the chain rule u — = 1 ) 

dx ) 


t \ d f u 

f(x) = m u . 

v ' dx l bL J 


The total energy of a body is 
me 2 

E = — — . 
bL 


Derivating Eq. (4) one gets 


. 2 A f uL 

dx I bL 


Let the right sides of Eq. (5) and Eq. (3) be equal, 


d 

2 1 ^ 

d 

u ^ 



= m n u — 


dx 

V. bLy 

0 dxl 

bL , 


Multiplying both sides of Eq. (6) by the factor dx, 

( I h ( u \ 

m n c d = m n u d . 

V bL J IbLj 


Integrating of Eq. (7), 


( 1 f f u | 

m n c 2 d — = m n ud . 

bL l bL J 
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The left hand side of Eq. (8) is the same as the total 
energy E. The integration by parts in the right hand side. 


E = m n u d 


\ ^ t \ l 

U U 2 U du 

— =m 0 + c - — — 

bL ) bL V c“ ) bL 


2 n , 2 

u , bL 


= mJ + c 


= m n + - 


bL bL 


2 2 2 \ „ 2 

u c u m n c 


l^bL bL bL J bL 
, N dE 

f (x) = — , ( 

dx 

when the left hand sides on Eq. (3) and Eq. (5) are equal. 
Simplifying the right sides of Eq. (2) , (3) , (5) one gets 


0 bL 3 ' 


Power = P, 
dE 


d [ 1 


dt l bL 


2 u du 1 

= m o c " ~T~7TT 
c dt bL 


= f -u 


In general the force is 


d { U 


F (t) = m 0 , 

dt l bL J 


d ( U 


dtl bL 


[ dU 2 1 UdU 1 

= m 0 bL 2 — + U-— - 

Vdt bL 3 c‘dtybL 3 

( dU U dU U dU 3 ) 1 

l dt c‘ dt c" dt ) bL 
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dt bL 

0 , c 2 , , dt , bL 3 


f dU ) 1 

= m o T -U 

l dt J bL 3 


= F • U 


or one can write 


dE 

F • U = . 

dt 

( 1 dt 

From (16) one gets — = 

U dR 


d f 1 


dRl bL 


U dU 1 


c 2 J l dR J bL 3 


= mJ ux I •- 


dux I 


- + uy J • 


duy J | 1 


dR ) bL' 


ux I, uy J are the component vectors of vector U, and are func 
tions of x, y. I, J are the unit vectors. 


Remembering dR = dx I + dy J, one gets 


f uxlduxl-l , x uyJduyJ-J , \ I 1 

F = (dx I + dy j) + — (dxl + dyj) — - 

^ I J y bL' 

dux duy^) 1 

= m 0 ux I + uyJ — - 

y dx dy JbL 


= m 0 (axl + ay j) — j 




dE 

F(R) = 

dR 


= grad E . 


( 20 ) 
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dU dU 

A = = U . 

dt dR 

The upgrowth of the kinetic energy in [Ri,R 2 ] : 


F • dR = nuc' 


!(i 

l 


du 

bL 3 


bL 


= e 2 -e 1 


when U-dU = udu; U-U = u 2 . 


( 21 ) 


( 22 ) 


In the last two files (The Quantization of Mathematics 1 and 
2) is used the polar coordinate system r,9 , and its conformal map- 
ping to the rectangular plane. The calculations typically go: 


^ dr du 

r,0: u = — : a = — 
d0 d0 


a = -u~'^(d 2 - u 2 y 7 

(23) 

j-u-\/(d 2 -u") 7 du 

(24) 

•J(d s -uY _ Q 
9 

(25) 

rl = jVd 2 -( 90) 2/9 d0 

(26) 

r2 = J— Vd 2 -(90) 2/9 d0 

(27) 


Simplify the right hand side of (26) giving (28) and 
simplify the right hand side of (27) giving (29) . 


(28) , (29) when substituted d = 5 can be plotted separately in some 
interval and in some scale, and they can be summed: ((28) + (29)) 

and then plotted. The sum of (28), (29) is analytically = 0, but with 
DERIVE ' s computing technique the program plots very interesting re- 
sults. DERIVE ' s version 2.51 does not do the same. 

The calculations, and especially plottings, take long time; I have 
plotted with my 486/33 two or three days. 


Why does the program plot the sum of (28), (29)? I think it makes use 
of roundings in calculations, which causes quantum effects in plot- 
ting. (28) and (29) are not identically the same, with reverse sign, 
and roundings are different in both. But together they form logical 
pictures in plottings, which I think the nature uses in the so cal- 
led quantum level. 


[ Erkki Ahonen: Kappaleen liikkeen kuvaaminen x,t-tasolla, 
Halikko 1983, ISBN 951-99465-7-8 ] 
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The first two examples : 

1: "THE STRAIGHT LINE" 

2 : "a-aO (u/uO) A pV (d A 2± (u/uO) A k) A n=0" 

3 : "a=d A 2x/dt A 2=du/dt=udu/dx : u=dx/dt" 

4 : "a0=0 ; a=0 ; u=uO " 

5 : "uO=Integration constant" 

6: "uOt^x-xO" 

7 : "xO=Integration constant" 

8: "u0=x0=0.5" 

9: l/2*t=x-l/2 

10: t=2*x-l 

11 : "plot straight line #10" 

12 : "u0t=x-x0" 

13: "uO = -0.5; xO = -0.5" 

14: -l/2*t=x+l/2 

15: t=-2*x-l 

16: "plot #16" 



Bitn tption 

Crtsi r:t S: Stilt x;2 y:2 hrive 26-plot 


Straight Line and Parabola 

1 : "THE PARABOLA" 

2 : "a-aO (u/uO) A pV (d A 2± (u/uO) A k) A n=0" 

3: "a=d A 2x/dt A 2 = du/dt = u du/dx: u =dx/dt" 

4: "p = n = 0; a = aO; a = du/dt" 

5: "du = aO dt: u = aO t" 

6: "dx = aO t dt: aO = 1" 

7: INT (t, t) 

8: t A 2/2=x 

9: t=-SQRT (2) *SQRT (x) 

10 : t=SQRT (2) *SQRT (x) 

11: "plot #9, #10" 

12: "p = n = 0; a = aO; a = du/dx" 

13: "u du = aO dx: u A 2/2 = aO x : aO = 1" 

14: u A 2/2=x 

15 : u=-SQRT (2 ) *SQRT (x) 

16 : u=SQRT (2) *SQRT (x) 

17: "dx/±V (2x) = dt" 

18: INT (1/SQRT (2*x) ,x) 

19: SQRT (2*x) =t 

20: -SQRT (2*x) =t 

21: "The equation of the parabola:" 

22: "aO (t - tO) = ±V(2 aO (x - xO) ) " 

23: "aO = 2;t0 = xO = 1/2" 

24: 2* (t-1/2) =SQRT (4* (x-1/2) ) 

25: t= (SQRT (2) *SQRT (2*x-l) +1) /2 

26: 2* (t-1/2) =-SQRT (4* (x-1/2) ) 

27: t= (1-SQRT (2) *SQRT (2*x-l) ) /2 

28: "plot #25, #27" 

Will be continued in the next issue! 

See a plot of Erkki ' s "Quantization" . 
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Trigonometry for the Classroom 

Josef Bohm, Wiirmla 

Yes. I know that the next contribution is easily to be done with any programming language. But. when 
I'm working for school or for private purposes I don't want to look for a while in my toolbox to find 
the fitting instrument. So, when I'm working with DERIVE, why don't use DERIVE for calculating 
triangles, too. 

On the other hand I think that this application can be easily done with students. I give the starting 
points and let them continue the work. We can train simple CASE-decisions, abstractions and at last 
the use of a selfmade tool. 

(sides si, s2, s3 with opposite angles al, a2 , a3 (in degrees) 


#1: Calculation of Triangles - Drei ecksberechnungen 

#2: the 3rd angle - der Erganzungswi nkel 

#3: R_A(al , a2) := 180° - al - a2 

#4: Triangle's inequation - Drei ecksungl ei chung 

TR_INEQU(sl , s2 , s3) := 

If si + s2 < s3 v s2 + s3 < si v si + s3 < s2 
#5: 1 

0 

#6: The Cosine Rule - Der Kosinussatz 



si- SIN(a2) 

#10: SS_0PPA2(sl, al, a2) := 

SIN(al) 



SA_0PPS2(sl, s2 , al, k) := 


1 


•(1 - 2 • k) + 180 -k 
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the sides sl,s2,s3 are known - die drei Seiten sl,s2,s3 sind gegeben: 

SSS(sl, s2 , s3) := 

If TR_INEQU(sl , s2 , s3) = 1 
"no triangle !" 

["Angle opposite SI: " , CA_OPPSl(sl , s2, s 3) ; "Angle opposite S2: ", 


CA_0PPSl(s2, s3, si); "Angle opposite S3: ", CA_0PPSl(s3, s2, si)] 
#14: The sides si, s2 and the angle a3 between them are known: 

#15: zwei Seiten SI, S2 und der ei ngeschl ossene Winkel A3 sind gegeben: 


Third side: 


CS_0PPA3(sl, s2 , a3) 


#16: SAS(sl , a3 , s2) := Angle opposite SI: CA_OPPSl(sl, s2, CS_0PPA3(sl, s2, 

. Angle opposite S2: CA_0PPSl(s2, si, CS_0PPA3(sl, s2, 


The side s3 and the adjacent angles al, a3 are known: 

eine Seite S3 und die beiden anliegenden Winkel Al, A2 sind gegeben: 


Third angle: 


R_A(al , a2) 


ASA(al , s3, a2) := 


Side opposite Al: SS_OPPA2(s3, R_A(al, a2) , al) 

Side opposite A2 : SS_OPPA2(s3, R_A(al, a2) , a2) 


zwei Seiten SI und S2 und der gegeniiberl i egende Winkel Al sind gegeben: 
hier konnen eine, keine oder zwei Losungen auftreten! 
two sides si and s2 and one opposite angle al are known: 
one, none or two solutions! 

zuerst eine Hi 1 fsfunkti on SSAH, dann die Fall unterschei dung : arbeite mit 
SSA(sl , s2 , al) : 


at first an auxiliary function SSAH, but work with SSA(...), CASE- 
deci si on : 
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Solution: Triangle 
Angle opposite S2: 

#26: SSAH(sl, s2, al, k) := 

Third angle: 

Third side: 
k + 1 

SA_0PPS2(sl, s2 , al, k) 

R_A(al , SA_0PPS2(sl , s2, al, k)°) 

1 ° 

SS_0PPA2(sl, al, R_A(al , SA_0PPS2(sl, s2, al, k)°)) . 

SSA(sl, s2, al) := 

If si > s2 

SSAH(sl, s2 , al, 0) 

#27: If s2 • SIN(al) / si > 1 

"no tri angl e ! " 

VECTOR(SSAH(sl , s2, al, k) , k, 0, 1) 

#28: Some examples: 

#29: a = 102, b = 61, c = 109 

#30: SSSC102, 61, 109) 

Angle opposite SI: 66.99036 

#31: Angle opposite S2: 33.39848 

. Angle opposite S3: 79.61111 . 

#32: a=100 , b = 160, (x = 33° 

#33: SSAC100, 160, 33°) 

T Solution: Triangle 1 1 [ Solution: Triangle 2 

Angle opposite S2: 60.62436 Angle opposite S2: 119.3756 

#34: , 

Third angle: 86.37563 Third angle: 27.62436 

Third side: 183.2406 . . Third side: 85.13396 

#35: b = 85.3, c = 23.4, <x = 58.6° 

#36: SASC85.3, 58.6°, 23.4) 

T Third side: 75.78758 ' 

106.1197 


#37: 


Angle opposite SI: 
Angle opposite S2: 


15.28023 
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#38: 

a=100 , b 

= 160, <x 

#39: 

SSA(100, 

160, 

53°) 

#40: 

no triangle ! 


#41: 

a=160 , b: 

=100, 

(X = 

#42: 

SSA(160 , 

100, 

33°) 


#43: 


Solution: Triangle 
Angle opposite S2: 
Third angle: 
Third side: 


1 

19.9012 
127.0987 
234.312 . 


Pentagon: 

AB = 30, BC = 52, AE = 30, 

AD = 65, BD = 60, EC = 68, ZEAB = 106,1° 

Find ED = d = ? and CD = c = ? 


D 



#44: Calculation of the pentagon - FLinfecksberechnung 

#45: [AB := 30, BC := 52, AE := 30, AD := 65, BD := 60, EC := 68, EAB := 106.1°] 


#46: EB = ? , AEB = ? 

#47: SAS(AE , EAB, AB) 


#48: 


#49: 


Third side: 
Angle opposite SI: 

. Angle opposite S2: 
[EB := 47.9496, AEB := 


47.94961 ' 
36.95 
36.95 . 

36.95°] 


#50: BEC = ? 

#51: SSSCBC, EB, EC) 


# 52 : 


Angle opposite SI 
Angle opposite S2 
. Angle opposite S3 


49.68443 ' 
44.67621 
85.63935 . 
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#53: BEC := 49.6844° 


#54: DAB = ?, EAD = ? 

#55: SSS(BD, AB, AD) 


#56: 


Angle opposite SI 
Angle opposite S2 
. Angle opposite S3 


66.98167 ' 
27.39935 
85.61897 . 


#57: EAD := EAB - 66.9816° 


#58: ED = d = ? 

#59: SAS(AE , EAD, AD) 


#60: 


Third side: 
Angle opposite SI: 
. Angle opposite S2: 


45.81712 ' 
24.40063 
116.4809 . 


#61: d = 45.8171 


#62: 

#63: 


DEC := 116.4809° - AEB - BEC 
SAS(45 . 8171 , DEC, EC) 


#64: 


Third side: 
Angle opposite SI: 
. Angle opposite S2: 


36.31207 ' 
38.89904 
111.2544 . 


#65: c = 36.3121 


I find it very useful to doublecheck the results by means of a plot in the 2D-plot window re- 
producing the sketch procedure how if it were done by hand or by using Dynamic Procedure. 
There are a lot of additonal advantages by transfering the abstract calculation to working with 
geometric objects and by connecting trigonometry with analytic geometry. I show the check 
performed with DERIVE and with my favourite Dynamic Geometry program GEOGEBRA 
which can be downloaded for free (for many languages) from www.qeoqebra.at . (2006) 
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#1: 

CaseMode := Sensitive 


#2: 

InputMode := Word 


#3: 

i 1 

> 

1 f 

i 1 

O 

0 

CO 

1 f 

1 — 1 
UJ 

o 

0]] 

#4: 

2 2 2 
kl := x + y = 65 


#5: 

r 2 2 2 

L kl := x + y = 65 , 

k2 := 


SOLUTIONSCkl a k2 , [x, y]) = 

' 305 5-720615 ' 

D := , 

.12 12 

[ 2 2 2 
L k3 := x + y = 30 , hi := y = TAN(106 

S0LUTI0NS(k3 a hi, [x, y]) 

T 8.319439599 -28.82337462 


2 2 21 
50) + y = 60 J 

305 5-720615 

12 12 

305 5-720615 

12 12 


.l°)-x. 


L -8.319439599 28.82337462 J 

#11: E := [-8.319439599, 28.82337462] 

2 2 2 2 
#12: k4 := (x - E ) + (y - E ) = 68 , k5 := (x - 30) 

1 2 

#13: S0LUTI0NS(k4 a k5 , [x, y]) 

' 58.00792189 43.81274142 ' 

#14: 

. -4.32736149 -39.05934271 . 

#15: C := [58.00792189, 43.81274142] 

#16: [A, B, C, D, E, A] 

#17: [ | E - D | , | C - D | ] = [45.81709477, 36.31215067] 


2 2 
y = 52 


See the simultaneous DERIVE plot on the next page! 
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The DERIVE plot: 



The English Geogebra construction (ci and di are the measured lengths): 




13 Free objects 

) 9 B = (30, 0) 

1 9 D = (25.38, 59.82) 

C3 Dependent objects 

( 5 * A = (0,0) 

9 C = (58.01,43.81) 

9 E = (-8.32, 28.82) 

•9 P = 2450.87 
-9 a: 28.82x + 8.32y= 0 
1 9 a i = 30 
1 9 b = 52 

J c: x 2 +y 2 = 4225 
d c 1 = 36.34 
v d: (x - 30) 2 + y 2 = 3600 
[ » d, = 45.79 

j 9 e, = 30 

f O g = 45.79 

| O h = 36.34 
: 9 a = 106.1° 
li Auxiliary objects 

9 e: (x - 30 ) 2 + y 2 = 2704 
1 9 f: (x + 8.32 ) 2 + (y - 28.82) 2 = 462 


• 

• 

ABC 

* 

2 

2 

V 
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Extended Fourier Series in DERIVE 

De Siquera & Domingos, Petropolis, Brasil 

1 - ABSTRACT 

Nowadays, the computers are very important to the progress of science and technology. Very 
complex algebraic calculations can be done by programs. We illustrate it with High Energy Physics 
and Engineering calculations. The algebraic calculations have showed to be a powerful tool for teach- 
ing Mathematics in the graduation courses also. Exhaustive calculations needed in the series expansion 
can be done quickly using algebraic computation. In this work we used DERIVE to implement the 
Fourier series with three functions in the interval. However, we cannot let the physical sense to escape 
from the students. In order to achieve this we solve a Strength of Materials Problem. 

2 - INTRODUCTION 

This work has an objective implement and extend an utility in the DERIVE program and also 
solves a problem in engineering Strength of Materials. This utility file implements the Fourier series 
until three functions in the interval. Nowadays, exists the Fourier series with just one function in the 
interval, in the utility file INT_APPS.MTH. 

Some series has been successfully tested and this stimulates the students of the Engineering cour- 
se to solve more elaborated problems. 



Figure 1 . A periodic function v(a) 


The period of the function in figure (1) is 2 p. This periodic function could be written as terms of sines 
and cosines as being (see [1]) 


1 ■A nnx A, . nnx 

y(x) = - a Q + 2^ a R cos + 2 b n sin 

^ n= 1 P n= 1 P 


(1) 


were a 0 , a„ and b„ are the coefficients of Euler or Euler-Fourier formulas and are: 
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( 2 ) 

( 3 ) 

( 4 ) 


where p is the semi-period. 

Substituting the equations (2), (3) and (4) in (1), we would obtain the Fourier series. 

There exists one theorem [2] that would not be showed here, that affirms, that for developing the 
Fourier series, the equations (2), (3) and (4) can be splitted in intervals [xi,x 2 ], |x 2 ,x 3 l, ...., [x„,x„+i], 
where each interval counts each function. 

In this work we will limit ourselves to just three functions in the interval [a,d\ to development of the 
Fourier series. These generalized three functions are showed in figure (2) 



Figure 2. A periodic function with 
three functions in the interval [a,d] 


Then the equations (2), (3) and (4) will be reduced to: 


j f b c d 

a Q = — I | (x) dx + 1 y 2 (x) dx + 1 y 3 (x) dx 


i f b c d \ 

lr . . nnx , f , . nnx , r , . nnx , 

a n = — \ jjtxjcos dx+ y,(x)cos dx+ y 3 (x)cos dx (7) 

pyi P b P c P ) 


b n = — [ f y 1 (x) sin dx + I* y , (x) sin dx + f >’ 3 (x) sin dx 

P\{ P t P t P ) 


( 8 ) 


with p showed in equation (5). 

Substituting the equations (6), (7) and (8) in (1) will result in Fourier series of the functions 
showed in figure (2). 
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4 - THE IMPLEMENTATION IN DERIVE 


For the implementation in DERIVE we define two functions, as showed below: 

F0URIER_2 FUNCTION (x, a, b, c, n, yl , y2) = RHS Equation (1) (9) 

F0URIER_3 FUNCTION (x,a,b,c,d,n,yl,y2,y3) = RHS Equation (2) (10) 

where: 

FOURIER_NFUNCTION ( ) - Name of the function defined in DERIVE; 

x - Variable that we intend to integrate; 
a, b, c, d - Borders of the intervals; 
n - Number of desired terms; 
y 1, yl, y3 - Functions in interval. 

By this way the work of expanding will be extremely reduced. 


5 - APPLICATION 

One application that grips engineering students was the fast solution of a complex problem of 
Strength of Materials, using the implementation in this work. 


Consider the beam support the 
load showed in figure (3) 



Figure 3. The beam to solution of engineering problem 


In the interval [0,7t] the distribution of the loading is sinoidal, between n and 2n it is uniformly dis- 
tributed and in \2n3n\ the loading is sinoidal again. Then for figure (3) we will write: 


y(*H 


1000(1 + sin x)Nm 1 
1000 Nm 1 


1 1000 (1 + sin x) Nm 1 


for 0 < x < 7t 
for n<x<2n 
for 2k < x < 3k 


( 11 ) 


Only one function could be used for describing the loading. Then, using the Fourier series and the 
implementation made in this work we can write: 


FOURIER 3 FUNCTION (x, 0, K, 2k, 3k, n, 1000 (1 + sin (x),1000,1000(l + sin(x)) 
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Figure (4) shows the developing of the Fourier series with n = 1 and n = 6. 



Figure 4. The loading developed in Fourier series 


The problem of Strength of Materials resumes in determining the bending moment diagram and 
the deflection of the beam. The solution for this problem would consist in solving the differential 
equation [3]. 

E I -j-j- — -W (x) (12) 

ax 

with E = Modulus of Elasticity; 

I = Moment of Inertia; 

W(x) = Loading. 

Equation (12) integrated two times gives us the bending moment: 


■ = —M (x) - - 1 dx | W (x) dx + c,x + c 2 


The boundaries conditions are: M(x = 0) = 0 and M(x = 3tt) = 0. 

Then we have c\ and C 2 - To solve equation (13) we define in DERIVE the function: 

MOMENT (x, cl, c2,W(x) ) = RHS Equation (13) (14) 

We will rapidly find Ci and c 2 equalling the above result to zero. 


The bending diagram is showed in figure (5). 

Figure 5. The bending moment diagram 

The maximum can be found by the graphic and 
your value is M = 14242 Nrn. 

(You can also differentiate and find the extremal 
value.) 



Integrating equation (13) two times we will find the function y(x) or the deflection of the beam. For 
this puipose we will define a new function in DERIVE: 
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Integrating equation (13) two times we will find the function y(x) or the deflection of the beam. For 
this purpose we will define a new function in DERIVE: 

- / dx (/ m dx + c3) + c4 

FUNCTION(x, c3 , c4, m) := 

E • I 


where C 3 and c 4 are constants that are fond by the boundaries conditions: 

FUNCTION^ = 0) = 0 and FUNCTIONS = 3m) = 0. 

(You will find the calculation at the end of the contribution. Editor) 

The deflection diagram is showed in Figure 6 with E = 210 GPa and I = 2.98 x 10 4 m 4 . 
(H-Beam W 250 x 167) 



Figure 6. The deflection of the beam 


The own weight of the beam was despized. We can find the maximum of y(x) by the graphic: 

Y(x)max = 0.0217 m. 

6 - CONCLUSION 

The main objective of this work was to extend the utilities for the DERIVE program, the calculation of 
the Fourier series up until three functions. It has been also shown one application on Strength of Mate- 
rials that would be almost impossible without the help of the DERIVE program. 

Some problems of Heat Transfer, Electricita, Waves would also use the Fourier series and they could 
be easily solved with the implementation showed in this work. 
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1 ( b c ' 

#1: F0URIER_2FUNCTI0N(x, a, b, c, n, yl, y2) := J yl dx + J y2 dx + 

c - a L a b , 


2 n ( 2-k-n-x 'j ‘ 

I COS • COS 

c - a k=l l c - a J . 


( b 

'j ( 2-k-n-x 'l 

COS • 

, J l c - a > 

. a 


f ( 2-k-n-x T 
yl dx + COS -y2 


' ( b c 

2 n ( 2-k-n-x f 2-k-n-x ~) 

dx + I SIN • SIN -yl dx + 

c - a k=l Lc-aJJ Lc-aJ 

, v a b 


C 2-k-n-x 'j 
SIN • y2 dx 


1 ( b c 

#2: FOURIER_3FUNCTION(x, a, b, c, d, n, yl, y2, y3) := J yl dx + J y2 

d - a L a b 


d ^ 2 n ( 2-k-7T-x ^ f ( 2-k-n-x ' 

dx + J y3 dx + Z COS • COS -yl dx + 

c J d - a k=l Ld-ayJ L d - a , 


' 2 • k-7T-x 'j ( 2-k-n-x \ 2 n 

COS -y2 dx + COS -y3 dx + I 

.. d-aj Ld-aJ d-a k=l 


' b c d 

' 2 • k-7T-x 'j f f 2-k-n-x ) f 2-k-n-x ) 

SIN • SIN -yl dx + SIN -y2 dx + 

d-aJJ ^ d - a J Ld-aJ 

^ a b c 


' 2 • k-7T'X 'j 

SIN • y3 dx 

. d-a 


FOURIER_3 FUNCTION (x, 0, n, 2-n, 3-n, 1, 1000- (1 + SIN(x)), 1000, 1000- (1 + 
SIN(x))) 


1200 -COS 


1000 ■ (3-n + 4) 


1200 -COS 


1000 ■ (3-n + 4) 


IF x > 0 a x < 3-n , 
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# 6 : 


#7: 


FOURI ER_3 FUNCTION (x, 0, n, 2-n, 3 -n, 2, 1000- (1 + SIN(x)), 1000, 1000- (1 
SIN(x))) 


6000 


f 4 ' x 1 f 2 -x I 

•COSl I 1200 -COSl I 


7 -77 


1000- (3-77 + 4) 


3-77 


#8: FOURIER_3FUNCTION(x , 0, 77, 2 -77, 3 -77, 6, 1000- (1 + SIN(x)), 1000, 1000- (1 

SIN(x))) 


6000 -COS 


#9: 1600-COS(4-x) 


10-x 


l 3 J 


1200 -COS 


8-x 


3 


9-n 


6000 -COS 


4-x 


L 3 


91-77 


1200 -COS 


2-x 


l 3 


11-n 


1000- (3-77 + 4) 


8000 • COS(2 • x) 


9-77 


7-77 


3-77 


#10: 


IF 


6000 -COS 


X > 0 A X < 3-77, 


8000-COS(2-x) 


1600-COS(4-x) 

9-77 
( 4-x 

6000 -COS 


10-x 


1200 -COS 


m 


l 3 


91 -TJ 
f 2-x 


1200 -COS 


V. 3 


11-77 


1000- (3-77 + 4) 


9-77 


7 -77 


3-77 


#11: FOURIER_3FUNCTION(x , 0, tt, 2-tt, 3-77, 12, 1000- (1 + SIN(x)), 1000, 1000- (1 + 

SIN(x))) 



Set w : =# 8 and let's proceed with this approximation for the given function 
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#12: w := FOURIER_3 FUNCTION (x, 0, n, 2-tt, 3 -n, 6, 1000- (1 + SIN(x)), 1000, 1000- (1 + 

SIN(x))) 

#13: MOMENT(x, cl, c2, w ) := J J w dx dx + cl-x + c2 


540 -COS 


#14: 100 -COS (4 -x) 


675 -COS 


2000-COS(2-x) 


3375 -COS 


2700 -COS 


2 2 
3-77-(500-x + cl-x + c2) + 2000-x 


#15: x = 0, M = 0 


#16: c2 


#17: c2 


3368935 


1716 -7T 


3368935 


1716 -7T 


#18: x = 3 tt, M = 0, c2 = 3368935 / (1716 tt)) 

#19: 3 - 7T - (1500 - 7T + cl + 2000) 

#20: cl = - 500- (3-77 + 4) 

#21: substitute for cl and c2: 


540 -COS 


#22: 100-COS(4-x) 


675 -COS 


2000-COS(2-x) 


3375 -COS 


2700 -COS 


f 2 3368935 

3 -7T- 500 • x + (- 500 - (3 • 7T + 4))-x + 

„ 1716-tt 


+ 2000-x 


Substitute for c i and C 2 to obtain m(.r)-expression #22, then set m: =#22 and plot the bending moment 
diagram (Figure 5). You can tty to find the extremal value. 


#23: m := #22 


#24: NSOLUTIONS — m = 0, x 

Idx 


#25: [4.712388980] 

#26: SUBST (m , x, 4.71238898) 

4 

#27: - 1.424206935-10 
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#28: FUNCTIONS, c3, c4, m) := 


J (J m dx + c3) dx + c4 


' 10 -x 
243 -COS 

#29: 25-COS(4-x) L 3 


6075 -COS 


500 • COS(2 • x) 


36-77-e-i 


30375 -COS 


455 -77- e • i 


6075 -COS 


2816-77-e-i 


9-77-e-i 


112-77-e-i 77-e-i 

2 3 4 3 2 

2574000-77 -x - 3432 • 77- (125 • x - 1000-x + 3-c3-x - 3-c4) - 5-x -(114400- 

10296-77-e-i 


x + 2021361) 


#30: Substitute for e = 210, i = 2.98 10 A 4 and the boundaries conditions: (— ^Figure 6) 

c4 350512647 

#31: 

6258000 381821440000 -n 


#32: c4 


1051537941 


183040 -7T 


3 2 

7T- (1287000 -7T + 1716000 -n - 1144 -c3 - 3368935) 

2386384000 
3 2 

5 • (257400 - tt + 343200 • n - 673787) 


#34: c3 


' 10 -x ) f 8-x 

81- COS 81-COS 

#35: C0S(4-x) l 3 J L 3 


C0S(2 • x) 


9011520-/T 949130000-tt 234967040-tt 112644-/T 


81- COS 


- cos ( i f) 


1869056-tt 


83440-77 


4 3 2 3 3 ~ 

1853280000-77 -x + 2471040000-77 -x - 411840000-77 -x + 2400 • 77- x- (28600 -x ~ 

103091788800- 


2 4 2 

228800 -x - 2021361) + 91520000 -x + 1617088800 -x - 9463841469 


00-77 
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#36: More examples: 

#37: Example 1: Two functions in the interval 

38: Original function: 

#39: [IF(x < -tt v x > tt, 0), IF(x > -tt a x < 0, 0), IF(x > 0 a x < tt, SIN(x))] 

#40: F0URIER_2FUNCTI0N(x , -tt, 0, tt, 3, 0, SIN(x)) 

2 • C0S(2 • x) SIN(x) 1 

#41: + + 

3 -TT 2 TT 



#42: Example 2: Three functions in the interval 


#43: 


' 

r 2 -tt 'i 


f 2 -tt 4-tt "1 

S 

LL 
1— 1 

x > 0 a x < , 2 

U_ 
1— 1 

X > A X < , 1 

, IF 


l 3 J 


L 3 3 J 

\ 


4-tt 

x > a x < 


2 -TT, 0 


#44: F0URIER_3 FUNCTION 


2-tt 4-tt 

x, 0, , , 2 -tt, 3, 2, 1, 0 


#45: 


3 • SIN(2 • x) 3 • SIN(x) 
+ + 1 


2-tt 


TT 


#46: F0URIER_3 FUNCTION 


2-tt 4-tt 

x, 0, , , 2 -tt, 8, 2, 1, 0 


#47: 


3 • SIN(8 • x) 3 • SIN(7 • x) 3-SIN(5-x) 3-SIN(4-x) 3-SIN(2-x) 

+ + + + + 


8-tt 

3 • SIN(x) 


7 -TT 


5 -TT 


4-tt 


2-tt 


+ 1 



TT 
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It is no problem to transfer the functions to the Tl-device. Don't include x into the variable list. 
It doesn't work. I show the Fourier-series for two pieces within the period and perform the 
major part of the calculation of the beam. Josef 

four_2 (a , b , c , n , y 1 , y2_) 

Func 

1/ (c-a)*(/(y1_,x,a,b)+/ (y2_, x , b , c) )+2 / (c-a) *I(cos (2*k*7r*x/ (c-a) ) * 

(/ (cos (2*k*7t*x/ (c-a))*y1_,x,a,b)+/ (cos (2*k*7t*x/ (c-a) ) * 
y2_, x , b , c) ) , k, 1 , n) + 21 (c-a) *1 (si n (2*k *k*x/ (c-a) ) * 

(/(sin (2*k*7i:*x/ (c-a))*y1_,x,a,b)+/(sin (2*k*7i*x/ (c-a) ) * 
y2_, x , b , c) ) , k, 1 , n) 

EndFunc 


tTi-TmoY ftp FH ’ r T fe T Ffi-r Y "i 

F f— |Fllgebra|CalG|ot,heriPrgnIo|ciean Up| | 


four-_2C - n , Q , ji , 3 , 0 , sintx)) 

■2-cos(2-x) sintx) . 1 
3 ■ ji 2 ji 

-2-cos(2-x) +i inM + i^ 2W Done 


3- ji 


'sinCx^Q < x and x<ji 
.. 0, else ' 


□one 


(0<x and x<n,sinCx>,G>-^yl<j<> 


RflP EH ACT 


FUNC 3 5 0 



TiT^Y fet y f:=t y fht V fs Y F b T Y i 
' 'iFUgebnalCalGlotherlPrgnlolciean Upl | 


m£ 


( undef , x < 0 and x > 3 ■ n 
f 10Q0 ■( 1 + sinCx)), x < ji 
l T 1G00, x < 2- ji 
| 1. 1000 ■( 1 + sinCx)), else’ 


else 


■» yl 

else |> 


□one 


■ four_3(0, ji, 2- ji, 3-ji, 6, 1000 (1 + sinCx))* 


-2QQ -1 SQQS ■ cos( 4 x) + 297Q - cos| ^L2L) + 


ans(l}->v2<x)l 


DNL05 F;h[i EHhCT 

FUNC E/30 


Pi=| 

fet Y n T FH T FET T FfiT Y F? T F :i 1 
Z o o r-i | T r a c e | R e g r a p h | M a t h | □ r a w | ▼ v 1 ! \ 

\ 

r\_f\ f 



DNL05 F;fi[i EHhCT FUNC 


ifl^oY FET Y f::t Y FHt T ™ y ft? y 1 

t— |fllgebra|Calc|0ther|PrgnI0|ciean Up| | 


■ y2(x> + i,j 

■ cos^-^j^-j -7-^18- cosj-^- ] + 5 (3 ■ ji + 4^ 


1 YX wdxtix 


- 500 (3 ■ ji + 4) ■ x^ - 3 ■ c ■ ji ■ x - 3 ■ d ■ Jijj 


<H,X,C> -x-d>l 


rri-TmoY ftp Y"F 3 *Hf FHT Y fe y fbt y i 

t— |fllgebra|CalG|0ther|PrgnI0|ciean Up| | 


Done 

1 nonentYo , d) = 0 | X = 0 

1716 ■ d ■ ji - 3368035 


1 soloe 


1716 ■ d • it - 3368835 
1716 ■ ji 


1716 ■ ji 
= 0, d] 
d = 


• = 0 


3368035 


1716-ji 


solue<ans<l>,d> 


n-7mnY FTP Y f3t y fit Y fe Y fbP y 1 
t f— |fllgebra|CalG|0t.her|PrgnI0|ciean Up| | 


■ soloe(3 (g + 500 (3 ■ ji + 4)) ■ n = 0 , g) 

g = -500 (3 ji + 4) 

■ nonentj “500 ( 3 ■ ji + 4) , Y^ie^if ) 

5- 30080 ■ gos( 4 ■ x) + 42768 gos| 10 3 * j + 11^ 


ans(l}-»v3<x>l 


DNL05 F;h[i EHhCT 

FUNC 13/30 
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Entdeckungen am Pascal 'schen Dreieck 

mit DERIVE 

Gunter Scheu, Pforzheim, Germany 

DERIVE is an outstanding tool to represent data for a conjecture 
quickly and systematically in a table or by a graphic. I want to 
show this on the Triangle of Pascal how to use DERIVE in this sense. 

Problems: Try to develop the term for one line in the T.o.P. 
Calculate the sum of the coefficients! 

What do you notice? 

What could be a rule for the sum? 

Proof your conjecture for the first 21 rows! 

And now for row 100! 

Produce a T.o.P with 11 rows! (or more!) 

Produce a T.o.P with rows' numbers and the rows ' s sums! 
Produce a T.o.P with a heading! 

Give the coefficients modulo 2, modulo 3, modulo k 

Consider the patterns! Try different values for k! 


#1: z(n) VECT0R(C0MB(n , k) , k, 0, n) 

#2: z(3) = [1, 3, 3, 1] 

#3: s(n) := I(z(n)) 

#4: s(3) = 8 

#5: VECT0R(s(n), n, 0, 20) 

#6: [1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, 16384, 32768, 

65536, 131072, 262144, 524288, 1048576] 

T 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 

18 

#7: Ll, 2, 2 , 2 , 2 , 2 , 2 , 2 , 2 , 2 , 2 ,2 ,2 ,2 ,2 ,2 ,2 ,2 ,2 


19 20] 

2 , 2 J 

r 100] 

#8: [s(100) , FACTOR(s(100))] = U267650600228229401496703205376 , 2 J 

#9: Di spl ayFormat:=Compressed 

#10: TABLE([z(n),s(n)],n,0, 11) 


0 [ 1 ] 1 

1 [ 1 , 1 ] 2 

2 [1,2,1] 4 

3 [1,3, 3,1] 8 

4 [1,4, 6, 4,1] 16 

#11: 5 [1,5,10,10,5,1] 32 

6 [1,6,15,20,15,6,1] 64 

7 [1,7,21,35,35,21,7,1] 128 

8 [1,8,28,56,70,56,28,8,1] 256 

9 [1,9,36,84,126,126,84,36,9,1] 512 

10 [1,10,45,120,210,252,210,120,45,10,1] 1024 


L 11 [1,11,55,165,330,462,462,330,165,55,11,1] 2048 J 
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#12: Z2b(n):rVECT0R(IF(M0D(C0MB(n,k) ,2)=0,“ 0 “),k,0,n) 


#13: TABLE(Z2b(n) , n, 0, 25) 


[ , ] 

[ , o , ] 


[ , 0 , 0 , 0 , ] 

[ , ,0,0, , ] 

[ , o , , 0 , , 0 , ] 


[ , 0 , 0 , 0 , 0 , 0 , 0 , 0 , ] 

[ , , 0 , 0 , 0 , 0 , 0 , 0 , , ] 


0 , 0 , 0 , 0 


23 [ 


, 0 , 0 , 0 , 0 , 0 , 0 , 0,0 


#16: Z2_b(n):=VECTOR(IF(MOD(COMB(n,k),2)=l, “1 

#17: TABLE(Z2_b(n) , n , 0 , 10) 


[ 1 , 1 ] 


[ 1 , 1 , 1,11 


J J J J 


[ 1 , 1 , , , 1 , 1 ] 


[ 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 ] 


8 [ 1 


J J J J J J J J 


9 [1,1 1,1] 


10 [ 1 , , 1 , 


J J J J J J J 
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#22: z3b(n):=VECT0R(IF(M0D(C0MB(n,k) ,3)*0, “ ,k,0,n) 

#23: TABLE(z3b(n),n,0,25) 


15 [ , % , 

16 [ , , % 

17 [ , , , 

18 [,%,%,% 

19 

20 

21 [ , % , % , , % , 
22 [,,%,,,% 

23 [ 

24 [ , % , % , , % , % , 

. 25 [ , , % , , , % , 


% , % , ] 


%,%,%, ] 




% , % , , % , % , ] 


% o/ o/ o/ o/ o/ o/ o/ 3 

, /o , /o , /o , /o , /o , /o , /o , J 


% 0 / 0 / 0 / 0 / 0 / 
I /o , /o , /o , /o , /o 






% 0 / 0 / 0 / 0 / 0 / 0 / 0 / 

, /o , , /o , /o , , /o , /o , , /o , /o , 

,%, , ,%, , ,%, , ,%, 


#25: pa(n , m):=VECTOR(IF(MOD(COMB(n , k) , m)#0 , “ & “),k,0,n) 

#26: pa_(n , m)"VECTOR(IF (MOD (COMB (n , k) , m)=0 , “ & “),k,0,n) 

#27 : bi ld(m , r)::TABLE(pa(l , m) , 1 , 0 , r) 

#28 : bi 1 dNEG(m, r):-TABLE(pa_(1 , m) , 1 , 0 , r) 

#29: [bi ld(3 , 10) , TABLE (z(n),n,0, 10) ] 


[ , ] 


# 30 : 5 


[ , & , & , ] 

[ , , & , , ] 


[ i & i&j ] 


7 

8 [ ] 

9 [ ] 


2 [1,2,1] 

3 [ 1 , 3 , 3 , 1 ] 

4 [ 1 , 4 , 6 , 4 , 1 ] 

5 [ 1 , 5 , 10 , 10 , 5 , 1 ] 

6 [ 1 , 6 , 15 , 20 , 15 , 6 , 1 ] 

7 [ 1 , 7 , 21 , 35 , 35 , 21 , 7 , 1 ] 

8 [ 1 , 8 , 28 , 56 , 70 , 56 , 28 , 8 , 1 ] 

9 [ 1 , 9 , 36 , 84 , 126 , 126 , 84 , 36 , 9 , 1 ] 

10 [ 1 , 10 , 45 , 120 , 210 , 252 , 210 , 120 , 45 , 10 , 1 ] 
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#35: bi 1 dNEG(4 , 27) 


24 [ & , , , , & , 

25 [ & , & , , , & , & 


#41: bildNEG(ll,27) 
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27 
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DERIVE's NEW FEATURES (last versions) 

1 . The new function RANDOM implements a pseudo-random number generator. Several 
functions using the random number generator are defined in the utility file MISC.MTH. 
(RANDOM_SIGN, RANDOMPOLY, RANDOMVECTOR, RANDOM_MATRIX). 

2. The functions MOD, FLOOR, GCD, and LCM formerly defined in the utility file MISC.MTH 
are now efficiently defined in DERIVE.EXE. FLOOR is now used to efficiently define the func- 
tions CEILING(m,n) (= [m/n]) and ROUND(m/n) in MISC.MTH. 

3. Number theory functions for finding the nth Fibonacci, Bernoulli, and Catalan numbers are 
now defined in MISC.MTH. 
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4. The function APPEND formerly defined in the utility file VECTOR. MTH is now efficiently 
defined in DERIVE.EXE. 

5. For convenience, the function FIT for finding the best least-squares fit to a parameterized 
expression now accepts two arguments: a label vector and a data matrix. 

6. The new function GOODNESS_OF_FIT defined in M ISC. MTH returns the goodness of the 
fit of the expression returned by the built-in function FIT. 

7. The Gamma function can now find a numerical approximation for complex as well as real 
arguments. 

8. The Riemann zeta function formerly defined in the utility file ZETA.MTH is now defined in 
DERIVE.EXE. In exact mode, ZETA(s) returns an exact closed form value if one is known. In 
approximate mode, it returns an approximate answer. Note that ZETA(s) is returned as the 
sum from k=1 to infinity of 1/k A s, if s>1. 

9. Functions for approximatig modified Bessel functions of the first and second kind are now 
defined in the utility file BESSEL. MTH. 

1 0. The utility files ODE1 .MTH and ODE2.MTH now include general purpose functions that 
automatically choose the appropriate method for solving a differential equation. Functions 
are provided for finding a specific solution when given an initial condition or boundary value, 
and for finding a general solution in terms of a symbolic constant c. 

1 1 . The utility file RECUREQN.MTH now has single-step functions for finding general and 
boundary value solutions to constant coefficient second-order linear difference equations. 

12. The new Transfer Load daTa command loads numeric data files produced by other ma- 
thematical software programs. 

13. Utility, demonstration, and numeric data files can now be loaded directly from DOS using 
a command line switch. 

14. The new Transfer Save C command saves expressions in a C language file in a form 
suitable for including in C programs. 

15. When a Plot command is issued from an Algebra window, you are now given the oppor- 
tunity to select the location for the new plot window. The choices are "Beside", "Under", or to 
"Overlay" the Algebra window. 

16. Powers of system or user-defined functions can be entered by putting the power immedi- 
ately after the function name. 

17. The new appendix "Line Editing Commands" in the DERIVE User Manual summarizes 
the line editing commands available when authoring an expression. Also the appendix "Addi- 
tional Resources" includes several new books based on DERIVE. 


In addition to the new features this version incorporates numerous internal improvements to 
DERIVE's factorizor, equation solver, and algebraic simplifier. 






